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Abstract 

In this paper we consider the effects on big bang nucleosynthesis (BBN) of 
the hadronic decay of a long-Uved massive particle. If high-energy hadrons are 
emitted near the BBN epoch (t ~ 10~^ - 10^ sec), they extraordinarily inter- 
convert the background nucleons each other even after the freeze-out time of 
the neutron to proton ratio. Then, produced light element abundances are 
changed, and that may result in a significant discrepancy between standard 
BBN and observations. Especially on the theoretical side, now we can obtain 
a lot of experimental data of hadrons and simulate the hadronic decay process 
executing the numerical code of the hadron fragmentation even in the high 
energy region where we have no experimental data. Using the computed light 
element abundances in the hadron-injection scenario, we derive a constraint 
on properties of such a particle by comparing our theoretical results with 
observations. 
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I. INTRODUCTION 



Big bang nucleosynthesis (BBN) is one of the most important tools to probe the early 
universe because it is sensitive to the condition of the universe from 10~^ sec to 10^ sec. 
Therefore, from the theoretical predictions we can indirectly check the history of the universe 
in such an early epoch and impose constraints on hypothetical particles by observational 
light element abundances. 

Now we have a lot of models of modern particle physics beyond the standard model, e.g., 
supergravity or superstring theory, which predict unstable massive particles with masses of 
O(IOO) GeV - 0(10) TeV, such as the gravitino, Polonyi field, and moduh. They have long 
lifetimes because their interactions are suppressed by inverse powers of the gravitational scale 
(oc 1/Mpi). These exotic particles may necessarily decay at about the BBN epoch (T < 0{1) 
MeV) if they have already existed in earlier stages. If the massive particles radiatively decay, 
the emitted high energy photons induce the electromagnetic cascade process. If the decay 
occurs after BBN starts t > 10^ sec, the light elements would be destroyed by the cascade 
photons and their abundances would be changed significantly. Comparing the theoretical 
prediction of light element abundances with the observations, we can impose constrains on 
the energy density, the mass, and the lifetime of the parent massive particle This 
subject was also studied in more details in the recent paper 

On the other hand, if the massive particles decay into quarks or gluons, near the BBN 
epoch 10"^ ^ t ^ 10^ sec. it is expected that the other important effects are induced. If once 
the high energy quarks or gluons are emitted, they quickly fragment into a lot of hadrons. 
Then, such high energy hadrons are injected into the electromagnetic thermal bath which 
is constituted by photons, electrons, and nucleons (protons and nucleons) at that time. At 
first, the high energy hadrons scatter off the background photons and electrons because they 
are more abundant than the background nucleons. Then, almost all kinetic energy of the 
hadrons are transfered into the thermal bath through the electromagnetic interaction. As a 
result, they are completely stopped and reach to the kinetic equilibrium. After that time, 
they scatter off the background p or n through the strong interaction, and they inter-convert 
the background p and n each other even after the usual freeze-out time of the neutron to 
proton ratio n/p oi the weak interaction. The effect extraordinarily tends to increase n/p. 
Therefore, the produced ^He would be increased in the hadron injection scenario compared 
to standard big-bang nucleosynthesis (SBBN). 

The pioneering investigation of this subject was done by Reno and Seckel (1988) and 
their treatments have been also applied to the other subjects 1^,0]. After their work was 
published, the experiments of the high energy physics have been widely developed. Now 
we can obtain a lot of experimental informations of the hadron fragmentation in the high 
energy region and also simulate the process even in the higher energies where we have no 
experimental data by executing the numerical code of the hadron fragmentation, e.g. JET- 
SET 7.4 Monte Carlo event generator In addition, we have more experimental data 
of the hadron-nucleon cross sections. Concerning BBN computations, it is recently needed 
that we perform a Monte Carlo simulation which includes the experimental errors of the 
reactions, and then we estimate the confidence levels (C.L.) by performing the Maximum 
Likelihood analysis and the fitting including both the theoretical and the observational 
errors. Performing the above procedures, we can compare each model in the various param- 
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eter sets. With these new developments in the theory, we set bounds to the hadronic decay 
of long-hved particles. [] 

This paper is organized as follows. In Sec. II we briefly review the current status of 
the observations and SBBN. In Sec. Ill we introduce the formulations and computations in 
the hadron injection scenario. In Sec. IV we compare the theoretical predictions with the 
observations. Section V is devoted to the summary and conclusions. 

II. CURRENT STATUS OF OBSERVATIONAL LIGHT ELEMENT 

ABUNDANCES AND SBBN 

A. Current status of observations 

In this section, we briefly summarize the current status of the observational light element 
abundances. The primordial D/H is measured in the high redshift QSO absorption systems. 
Recently a new deuterium data was obtained from observation of QSO HS 0105+1619 at z 
= 2.536 |[l^. The reported value of the deuterium abundance was relatively low, (D/H)"^* = 
(2.54 ± 0.23) X 10^^. Combined with the previous "low D" data |]TT|, the authors reported 
that the primordial abundance is 

lowD : (D/Uy'" = (3.0 ± 0.4) x 10^^ (1) 

We call this value "low D." On the other hand, Webb et al. obtained high deu- 
terium abundance in relatively low redshift absorption systems at z = 0.701 towards QSO 



PG1718+4807 |T2| 



highD : (D/H)"'" = (2.0 ± 0.5) x 10"^ (2) 



In these days, Kirkman et al. |TB[ also observed the clouds independently and obtained new 
spectra using HST. They claimed that the absorption was not deuterium although there were 
still some uncertainties. Here we think that it is premature to decide which component is 
correctly primordial and the possibility of "high D" have not been excluded yet. Therefore, 
we also consider the possibility of "high D" and include it in our analysis. 

The primordial value of ^He is inferred from the recombination lines from the low metal- 
licity extragalactic HII regions. The primordial value of ^He mass fraction Y is obtained to 
regress to the zero metallicity 0/H — for the observational data because it is produced 
with oxygen in stars. In these days Fields and Olive reanalyzed the data including the Hel 



absorption effect Then they obtain the observational V, 

Y°'' = 0.238 ± (0.002),^,^ ± (0.005),^,,*, (3) 

where the first error is the statistical uncertainty and the second error is the systematic one. 
We adopt this value as the observational value of Y. 



^For relatively longer lifetimes, there exists an another interesting process that the emitted high 
energy nucleons destroy the light elements which have already been produced 0]. 
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It is widely believed that the primordial abundance of ^Li/H is observed in the Pop II 
old halo stars whose temperature is high T^s ^ 6000 K and metallicity is low [Fe/H] < - 1.5. 
They have the "plateau" structure of ^Li/H as a function of the metallicity. We adopt the 



recent measurements by Bonifacio and Molaro 



log 



10 



('Li/H 



obs 



-9.76 ± (0.012),^,^ ± (0.05),,,^ ± (0.3) 



I add- 



(4) 



Here we have added the additional uncertainty for fear that the ^Li in halo stars might have 
been supplemented (by production in cosmic-ray interactions) or depleted (in stars) 



B. Current status of SBBN 

Here we show the current status of standard big-bang nucleosynthesis (SBBN). Within 
recent years, there was a great progress in the experiments of the low energy cross sections 
for 86 charged-particle reactions by the NACRE collaboration |T^. In the compilation, 



22 reactions are relevant to the primordial nucleosynthesis, and the old data were revised. 
In particular, 7 reactions of them are important for the most elementary processes up to 
mass-7 elements. On the other hand, recently Cyburt, Fields and Olive reanalyzed the 
NACRE data and properly derived the la uncertainty as a statistical meaning and the 
renormalization of the center value for each reaction [EO]. In addition, they also reanalyzed 



the four remaining reactions, using the existing data pT|-p3|] and the theoretical prediction 
(for one reaction) Their efforts are quite useful for the study of the Monte Carlo 

simulation in BBN, and it was shown that their treatment is consistent with the other 
earlier studies adopting the results of NACRE [PHEB . 



Carrying the Monte Carlo simulation into execution, we adopt the theoretical errors and 
the center values for 11 elementary nuclear reactions in Ref. For the error and the 

center value of neutron lifetime, we adopt the compilation of Particle Date Group p^, see 
Eq. (|To|). To systematically take account of the uncertainties, we perform the Maximum 
Likelihood analysis including both the observational and theoretical errors which are 
obtained in Monte Carlo simulation. Here we assume that the theoretical predictions of 
(D/H)*^, y*^, logio[('^Li/H)*''] obey the Gaussian probability distribution functions (p.d.f.'s) 
with the widths given by the la errors. Concerning the observational values, (D/H)"'"', 
and logio[(''Li/H)°'''^] are also assumed to obey the Gaussian p.d.f.'s. 

In Fig. we plot ^is a function of baryon to photon ratio, t] = UB/n^, where hb is the 
baryon number, and is the photon number. The solid line (dashed line) represents the 
case of low D (high D). From this figure, we find that SBBN agrees with the observation 
of ^He, D, and ^Li very well at more than 95 % C.L., and we obtain i] = 5.6lo;8 x 10~^° 



^These days, however, it was claimed that there is a significant Li-Fe trend in the low metallicity 
region [|l^. In addition, Ryan et al. ||l8| assumed that this trend is due to the cosmic ray interac- 
tions, and they inferred the primordial value is ^Li/H = (1.23^032) x 10"^'^. Because we can not 
make a judgment on the above discussions, for the moment we adopt the value in Eq. (^) with 
large uncertainties in this paper. 
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(r] = l.SZ'oi X 10"^°) for fow D (high D) at 95 % C.L. Using the relation Qsh^ = 3.63 x 
Wr]{To/2.725K), we obtain 



U.UZUO_o.o029 



U.UUUd_0.ooi8 



(5) 



(for lowD), 
(for highD), 

at 95 % C.L., where fi^ is baryon density parameter, h is normahzed Hubble parameter as 
-f^o = lOO/i km/sec/Mpc, and Tq is the present temperature [27|. Under these circumstances, 
we can check the non-standard scenario comparing the predictions of the BBN computations 
with observations. 



III. HADRONIC DECAY AND BBN 

In this section, we discuss the hadron-injection effects on the history of the universe 
near BBN epoch {t = 10^^ - 10'^ sec). Here we consider the case that the unstable massive 
particle "X" has some decay modes into quarks and gluons, and as a result it induces the 
late-time hadron injection. 



A. Time scale of the interactions 

If the quarks and gluons were emitted by the decay of the parent particle X whose mass 
is about O(IOO) GeV - 0(10) TeV, they immediately fragment into hadron jets and produce 
a lot of mesons and baryons (vr^, vr", K^, g, n,p, A", and so on). Then, the typical energy 
of the produced hadrons is about C(l)GeV - C(100)GeV, and they are injected into the 
electromagnetic thermal bath which is constituted by 7,e^, and nucleons. 

As we see later, if once such high energy hadrons are injected into the thermal bath in the 
beginning of the BBN epoch (i.e., at the temperature T ^ 0.09 MeV), their almost all kinetic 
energy is transfered into the thermal bath through the electromagnetic interactions except 
for neutral kaons. Then, such hadrons scatter off the background particles, and then they 
induce some non-standard effects on BBN. Especially, the emitted hadrons extraordinarily 
inter-convert the ambient protons and neutrons each other through the strong interaction 
even after the freeze-out time of the neutron to proton ratio n/p. For the relatively short 
lifetime {r^ ^ 10^^ sec - 10^ sec) in which we are interested, the above effect induces the 
significant change in the produced light elements. Concretely, protons which are more 
abundant than neutrons, are changed into neutrons through the hadron-proton collisions 
and the ratio n/p increases extremely. In this case, the late-time hadron injection scenario 
tends to increase ^He because it is the most sensitive to the freeze out value of n/p, 

The emitted hadrons do not scatter off the background nucleons directly. At first hadrons 
scatter off the background photons and electrons because they are much more abundant than 
background nucleons (about 10^° times larger). As we see later, for t < 200 sec, the emitted 
high energy hadrons are immediately thermalized through the electromagnetic scattering 
and they reach to the kinetic equilibrium before they interact with the ambient protons, 
neutrons and light elements. Then we use the thermal-averaged cross section (erf )^i^jy, for 
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the strong interaction process N+Hi N'-\ — between hadron Hi and the ambient nucleon 
A^, where denotes proton p or neutron n. The strong interaction rate is estimated by 



where n^v is the number density of the nucleon species N, t] is the baryon to photon ratio 
(= ^B/n^), riB denotes the baryon number density (= np + n„), and is the nucleon 
fraction (= ni^/uB)- Here, for the moment we adopt 40 mb as a typical value of the cross 
section for the strong interaction. This equation shows that every hadron whose lifetime 
is longer than (9(10~^) sec contributes to the inter-converting interaction between neutron 
and proton at the beginning of BBN. Hereafter we will consider only the long-lived mesons 
( TT^, K^, and Kl) and baryons (p, p, n, and n). ^heir lifetimes are given by 



r^± = (2.6033 ± 0.0005) x 10"^ sec, (7) 

Tk± = (1.2386 ± 0.0024) x 10"*^ sec, (8) 

r^o = (5.17 ± 0.04) X 10"® sec, (9) 

r„ = 886.7 ± 1.9 sec, (10) 

and proton is stable. 

Here we define the stopping time Xg^p of the high energy particle Hi in the thermal 
plasma as 




(11) 

where E denotes the energy, dE / dt denotes the energy loss rate in the thermal plasma and 
it depends on the each scattering process of particle Hi off the background particles. Eq 
is the initial energy, and i?th is the threshold energy of the process, ^o estimate whether 
particle Hi is stopped or not in the thermal plasma through the electromagnetic interaction 
until it scatter off the background baryons (n, p, and produced light elements), we computes 
the rate, 

-^stop = '^N-iN' ^ ''"stop' (12) 

as an indicator which roughly represents the number of the scattering during the stopping 
time T^op. If -R^op is much less than unity, the emitted high energy hadron Hi is completely 



^TT^,Kg, and have much shorter lifetimes and they have completely finished to decay because 
their lifetimes are r^o = (8.4 ± 0.6) x IQ-^^sec, r^o = 0.89 x 10"^°, and tj^o = 2.63 x 10"^° sec 
respectively. Therefore, they do not contribute to the interesting process in this situation. 

'^To roughly estimate the timescale until the particle is stopped, it would be usually adequate that 
we take -Eth to be equal to the mass of the particle Hi in the relativistic regime. 
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stopped and can not reach the background baryons with the high energy. On the other 
hand, if -R^op is greater than unity, the high energy hadron can not be stopped through the 
electromagnetic interaction and directly scatters off the background baryons. In addition, it 
might destroy the light elements which have already produced if the particle X decays after 
the cosmic time is t ~ 200 sec. 



B. Hadron stopping in the electromagnetic thermal plasma 

When the cosmic temperature T is higher than the electron mass mg, there are sufficient 
electrons and positrons in the universe. In this situation, it is expected that the emitted 
charged particles vr^, K"^, and p are quickly thermalized through the electromagnetic inter- 
action. In fact, the energy loss rate of the charged particle through the Coulomb scattering 
is given by 

for T ^ rrie in the relativistic regime, a is fine structure constant (~ 1/137). Then, the 
stopping time of the charged particle ( "ch" ) is estimated by 



r.t„.1.18xlO-Sec(^)(^)"\ (14) 



for T ^ me. Then, -Rg^p is much smaller than unity and we can regard that charged hadrons 
are completely stopped. 

As for neutron, we can see that it is also completely stopped for T ^ me- Although 
neutron is neutral, of course, it can scatter off the background electrons through the elec- 
tromagnetic interaction by the magnetic dipole moment. The energy loss rate through the 
Coulomb scattering is given by 

dE _ 15m^ 

dt rn^a^glT*' ^ ^ 

in the relativistic regime, where Qn is neutron magnetic moment (= —1.913) |2^, and m„ is 
neutron mass. The stopping time of neutron is 

r°„p ^2.34x10-1° sec (^). (16) 

Thus -Rstop of neutron is much smaller than unity, and it does not scatter off the background 
baryons before it stops for T ^ mg. 

On the other hand, if the temperature is much lower than electron mass (T < rrie), the 
situation is quite different because the number density of electrons becomes little. In this 



^Although the above estimations have been discussed only in the relativistic regime, the similar 
results are also obtained in the non-relativistic regime B. 
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case, the emitted mesons completely decay and disappear in the universe before they scatter 
off the background baryons because the lifetime is shorter than the timescale of the strong 
interaction (see Eq. (^). Thus we should not treat the injection of any mesons in such a late 
epoch. Because proton is stable, and neutron has a long lifetime compared to the typical 
timescale of the strong interaction in Eq. (^), we should worry about the thermalization of 
the emitted high-energy nucleons till quite later. 

For proton, the ionization loss is more effective to lose the relativistic energy for T ^ me- 
The ionization-loss rate is expressed by 

dt V P \ LUp J ^ ^ 

where Z denotes the charge {Z = 1 for proton), v is the velocity of the high energy proton, 
7 is the Lorentz factor, A is 0{1) constant, and ujp denotes the plasma frequency. 



2 4:7rane , . 

% = ) (18) 



where Ue represents the electron number density. We evaluate the stopping time of the 
proton to lose its relativistic energy. 



''"stop 



1.2xl0-»=ec.^e^(^)(f)-'. (19) 



where 7]iq is defined hj rj = rjiQ x 10^^'^, and the dimensionless parameter x = rrie/T. If we 
demand -Rstop ~ 1, we obtain T > 22keV which corresponds to cosmic time t ^ 3 x 10^ sec. 
Namely after t ~ 3 x 10^ sec, such a high energy proton can not be stopped in the thermal 
bath, and it is inevitable to scatter off the ambient baryons with the high energy. 

As well as the high energy proton, we estimate the case of the high energy neutron. The 
energy loss rate of the neutron through the Coulomb scattering for T ^ me is 



dE Sna'^g'^m, 
dt 



neE\ (20) 



n 



The stopping time to lose the relativistic energy is 

r,'^„p^l.68xl0-^ecxie^(^) . (21) 

Here if we require -R"top ^ 1, we find that the temperature should be greater than 95 keV for 
the neutron stopping which corresponds to the condition that cosmic time should be shorter 
than 150 sec. In this case, after t ~ 150 sec, the high energy neutron essentially inevitably 
scatters off the background baryons before it stops. Under these situations, at longest, after 
t ~ 150 sec the high energy nucleons necessarily scatter off the ambient baryons through the 
strong interaction, and we would also have to worry about the possibilities of the destruction 
of the light elements. This means that the scattering process after t ~ 150 sec is beyond the 
limits of validity in our treatment in this paper. For the problem, we will discuss it later. 

As for it is never stopped in the electromagnetic plasma because it does not interact 
with electrons and photons. Therefore, by using the energy dependent cross sections we will 
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treat the scattering off the ambient nucleons. To perform the computation, we should know 
the correct energy distribution of produced through the hadron fragmentation. 

On the other hand, for relatively longer lifetimes Tx ^ 10^ sec, there is another interesting 
effects on BBN. The emitted photons or charged leptons induce the electro-magnetic cascade 
showers and produce many soft photons, ^heir spectrum has a cutoff at E™^^ ^ ml/ (22T). 
If E^^^ exceeds the binding energies of the light elements, these photons destroy the light 
elements and change their abundances 0,^. In fact, at t ^ 10^(10^) sec, the energy of 
the photon spectrum which are produced by the decay of X, exceeds the deuterium C^He) 
binding energy B2 — 2.2 (i?4 ~ 20) MeV. However, because we are not interested in the 
photodissociation here, we only study the case of tx ~ 10^ sec. 

C. Hadron Jets and collider experiments 

As an example of the hadronic decay, if the gravitino V'^i is the parent particle X whose 
mass is mx = 0{1) TeV, it can have net hadronic decay modes, e.g., tp^j, ^qq (g: quark), 
with the branching ratio Bh- In this case, Bh can become ~ 0{a) at least even if the main 
decay mode is only — 77 (7 : photino), because of the electromagnetic coupling of the 
photon. As we quantitatively show later, about 1 hadrons are produced for B^ = 0.01 and 
for the energy per two hadron jets, 2Ejet ~ 2/3mx, if we assume that the mechanism of the 
hadron fragmentation is similar to the e^e~ collider experiments. In addition, the emitted 
high energy photon whose energy is about ~ mx/2 scatters off the background photon 7bg 
and can produce a quark-antiquark pair. [jThen, the center of mass energy is about ^/s ~ 2 
GeV and produces about 3 hadrons which could effectively contribute to the decay mode into 
hadrons as the branching ratio Bh ~ (9(0.01). Therefore, we should consider the hadronic 
decay modes at least as Bh = (9(0.01) in this case. On the other hand, if the decay mode 

~^ 99 {9- gluon, and g : gluino) is kinematically allowed, Bh may become close to one. 

For the other candidate of the parent particle, Polonyi field or moduli, which appears 
in supergravity or superstring theory and has a (9(l)TeV mass, would also have a hadronic 
decay mode {(f) ^ gg ). 

Fortunately, we can estimate the number and energy distribution of the produced hadrons 
by using the JETSET 7.4 Monte Carlo event generator [§]. This FORTRAN package com- 
putes the hadron fragmentation for the qq event (g: quark) in the e~^e~ annihilation and 
predicts the energy distribution of the products to agree with the e~^e~ collider experiments. 
In Fig. we plot the averaged charged-particle multiplicity {Nch) which represents the total 



^Even if the decay modes into hadrons are dominant (Bh ~ 0{1)), the almost all parts of the 
energy of the parent particle are transfered into photons and electrons because the hadrons decay 
after they completely transfer their relativistic energy into the thermal bath. In addition, it is 
expected that about 1 /3 parts of the produced hadrons are approximately vr'^ and they decay as 
Tj-O with a much shorter lifetime (r^o ~ 10^^^ sec). 

^Of com'se, there are some leptonic modes in the process, e.g. 7 + 7bg — > e"*" + e~. Thus, the net 
branching ratio into hadrons is about ~ 60% in this energy. 
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number of the charged hadrons emitted per e~^e annihilation and per two hadron jets as a 
function oi ^/s {= 2 Ejet). 0Recently LEPII experiments (ALEPH, DELPHI, L3 and OPAL) 
give us the useful data for = 130 - 183 GeV. Therefore, now a number of experimental 
data are available at least up to y/s ~ 183 GeV . The filled circle denotes the data points 
of e~^e~ collider experiments. From Fig. |^ we find that the predicted (A^^ch) excellently agrees 
with the experimental values. Thus, in this situation we use the JETSET 7.4 to infer the 
spectrum of the emitted hadrons extrapolating to the various higher energies. 

In Fig.0 we plot the spectrum of the produced mesons (tt"*" + tt", + , and A'^) 
as a function of the kinetic energy -Ekin- This is the case that the center of mass energy is 
y/s = 91.2 GeV which corresponds to the resonance. In similar fashion, in Fig. ^ we 
plot the spectrum of the produced baryons ((a) n + n, and (b) p + p). In Fig. |^ we plot the 
averaged number of the produced hadron per two hadron jets as a function of 2i?jet, which 
is obtained by summing up the energy distribution. From Fig. H, we find that almost all 
hadrons are composed of pions. 



D. Cross sections of hadron- nucleon scattering 

Because in this paper we are interested in the BBN epoch, i.e. T ^ 0{1) MeV, the 
temperature is much less than the typical mass of the emitted hadrons, e.g. rriHi = O(IOO) 
MeV - 0{1) GeV. As we discussed in section |IIIB| , as long as the temperature is relatively 



high enough (T > 95keV), the emitted high energy hadrons (tt^, AT^, p, and n) have com- 
pletely lost their relativistic energies through the electromagnetic interaction in the thermal 
plasma and are quickly thermalized except for neutral kaon K^. Then only the exothermic 
process is relevant for the hadron to scatter off the background baryons through the strong 
interaction because it has just a little kinetic energy of the order of the temperature T. Of 
course, such a low energy hadron can not destroy the background "^He. Concerning exother- 
mic reactions, it is well-known that the cross section a is nearly inversely proportional to 
the velocity v of the projectile particle in the low energy. Namely av almost does not have 
a V dependence and is nearly a constant for the beam energy. Except for A'^, therefore, we 
can use the threshold cross section instead of the thermal-averaged cross section. Here we 
adopt the results of the thermal-averaged cross section in Ref. ||^. 
The thermally averaged cross sections for tt^ are 

{av):l^ = 1.7 mb, (22) 
{cTv);Zn = l-5a(r) mb, (23) 

where Ch^ (T) is the Coulomb correction factor when the beam particle Hi is the charged one. 
Because the reaction p"*" + tt^ ^ n + ■ ■ ■ is enhanced due to the opposite-sign charge of the 



^Here {N^h) is defined as the value after both Ks and A'' had completely finished to decay, where 
their lifetimes are r^^o = 0.89 x 10"^" sec and t/^o = 2.63 x lO"^'^ sec respectively. As we showed 
in section [II A , we should not treat any particles with the shorter lifetime than ~ 10~^ sec in this 
situation. 
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initial state particles, we should correct the strong interaction rates by simply multiplying 
CHi{T) to that which are obtained by ignoring the Coulomb corrections. The Coulomb 
correction factor is generally estimated by 

g".'^)= l!'l^'S,T| . (24) 



where ii{T) = Q;y/ij/2T, a is the fine structure constant, and fii is the reduced mass of the 
hadron Hi and the nucleon. 

The thermally-averaged cross sections for K~ are 

{av)'^:^^ = 26 mb, (25) 

{av)^Zn = 34 mb, (26) 

{av)^Zn = 31Ck{T) mb, (27) 

{av)^Zp = U.5Ck{T) mb. (28) 

Here we ignore interaction because n + —>■ p + is the endothermic reaction which 
has Q = 2.8 MeV, and it is expected that the kinetic energy of is less than Q. 

As for neutral kaon K^, there are no adequate experimental data of the differential cross 
sections as a function of the beam energy to use in our current purpose. It is very serious 
for us because does not lose their relativistic energy and is never stopped in the thermal 
bath, and then we should know the differential cross sections in whole relevant energy range. 
For example, in Fig. ^ we find that the source distribution function of is spread in the 
wide energy range. At least we want to obtain the data of the cross sections for the typical 
i^°-beam energy, -Ebeam = 10 MeV - 1 TeV, where -Ebeam is the kinetic energy of K^^. In 
this situation, we should estimate the data table of the cross sections of the scattering 
by using the other existing experimental informations. 

Here we assume that scatters off the nucleon as a combination of 1/2 and 
1/2 K because in fact A'^ is nearly the linear combination of and K states that 
\Kl) ~ (lAfO) - \1C))/V2. gin addition, we assume that the strangeness of (K ) is 
similar to K"^ {K~) because A'" = ds {K^ = sd) contains s (s) {s: strange quark, and d: 
down quark). Of course, the above assumption is not wrong very much under the isospin 
SU(2) transformation for the doublet {u: up quark) because we cannot imagine that 
there exists a special coherent interference in the inelastic scattering. 

In this assumption, we would also have to worry about the effect of the Coulomb correc- 
tions because actually A'^A^ scatterings are not supposed to suffer from any electric charges. 
From Eq. (|^), however, we find that the Coulomb correction is less than 10% at most in 
both cases of the attractive force and the repulsive one as long as the kinetic energy of 
is more than (9(10) MeV. Therefore, we can ignore the Coulomb correction and the above 
assumption would be reasonable in this situation. 



^Of course, the CP violation effect does not change our rough estimates at all and is not important 
here. 
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Fortunately, we have good compilations of the experiments for the total cross section 



and the elastic cross section for K^p and K p [£7|- Thus, by averaging them we can 
estimate the total cr^Op and elastic cross section cr^^Op respectively. In Fig. |^, we plot the 

obtained total and elastic cross sections for K^p scattering. It is happy that the obtained 
total cross sections agree with the direct experimental data and the theoretical predictions 
marginally within a few tens percent although they were studied only in the low energy 



regions for -Ebeam ~ 0.7 GeV In addition, we have the experimental data of the inelastic 



scatterings, K^^p —>■ K^p and K^^p KgpTc^n [2£] which are also plotted in Fig. |^. Now 
we assume that the cross section of the inter-converting reaction K^^ + p —>■ n + ■ ■ ■ is 

obtained by apJin = l/2[a^°l^ - (o-^Op^^Op + ctrIp^kIp + (^kIp^k%pit+tv-)] because the final 
states of the inelastic scattering without K^^p — > K^p + ■ ■ ■ are KNtt, A^vr, or Svr, and it is 
approximately expected that either p or n appears in a closely even probability. QThen, we 



get the remaining cross section as ap-^p = cr^Op — c"; 



_ ^tot ^^L 



'p- 



About neutron- i^'^ scattering, we could have performed the similar treatments. However, 
compared to the cases of proton, we do not have adequate compilations for the neutron- i^'^ 
process. On the other hand, the data tell us that we can approximately regard them as the 
cross sections of the proton- ii" scattering within a few tens percent in the high beam energies 
{Ek ~ 100 MeV). The theoretical reason is that the strong interaction does not distinguish 
between proton and neutron in such a high energy. Under these circumstances, we assume 
that the cross section of K^j^n is as same as K^p with a few tens percent error. 

To perform the numerical computations including the i^^ -injection effects in BBN, it is 
useful to average the cross sections by the energy spectrum of K^. As we discussed in the 
previous subsection, we can use the JETSET 7.4 Monte Carlo generator and get the energy 
spectrum of emitted in wide range of the source energy. For example, we can see the 
spectrum of the produced for various energies in Fig. ^. Then, we get the averaged 

cross sections, (Tp4p and Op-^n, as the convolutions of the data of the cross sections with the 
energy spectrum of (Fig. |^). 

Concerning the emitted nucleons, we basically follow the Reno and Seckel's treatment 
that we regard the nucleon-antinucleon pair as a kind of a meson A^A^ p. Then, the A^A^ 
meson induces the inter-conversion A^ + A^A^ — > A^' + ■ ■ - . In Ref. we have the thermally- 
averaged cross sections, 

{av)f_^^ = 37 mb, (29) 
{av)f^^ = 28 mb, (30) 
{avYjUp = 28 mb, (31) 
{av)f_^p = 37 mb. (32) 

As we discussed in the previous sections, however, the late time emission of the high energy 



^OThe branching ratios are presented as n7r°(35%), p7r-(63.9%); S° A07(100%); S+ 

n7r"'"(48.3%), ^-^""(51. 6%); S~ — > n7r~(99.9%) [27|. We also ignore the multiple production process 
of baryons because the center of mass energy is too low for the process to dominate the other 
reactions. 
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nucleons would induce the destruction of light elements for T ^ 95keV. However, for the 
moment we treat the nucleons as if they are approximately thermalized. We will also discuss 
the modification on the result caused by the above simple assumption later. 



E. Formulation in hadron-injection scenario 

We formulate the time evolution equations in the late-time hadron-injection scenario 
here. As we have mentioned in the previous subsections, the hadron injection at the begin- 
ning of BBN enhances the inter-converting interactions between neutron and proton equally 
and the freeze out value of n/p is extremely increased. Then the time evolution equations 
for the number density of a nucleon N{= p,n) is represented by 

driN 



+ 3H{t)nN 



dt 



dt 



- BhVxUx {Kn^n' - Kn'-,n) , (33) 

- weak 

where H{t) is Hubble expansion rate, [dnx / dt]^^,^^^ denotes the contribution from the usual 
weak interaction rates as well as SBBN, is the branching ratio of the hadronic decay 
mode of X, nx is the number density of X, K^^n/ denotes the average number of the 
transition N ^ N' per one X decay. 

The average number of the transition N N' per one X decay is expressed by 

Kn^n'=T.^N^'R^^n', (34) 

where Hi runs the hadron species which are relevant to the nucleon inter-converting reac- 
tions, Xjet is the number of the hadron jet per one X decay, denotes the average number 
of the emitted hadron species Hi per one X decay. is presented in Fig. H as a function 
of 2£^jet, where i^jet is the energy of a hadron jet. Rj/^j.^, denotes the probability that a 
hadron species Hi induces the nucleon transition X X' and is represented by 

dH, _ ^N-.N' (OK) 
dec abs 

where = t^^ is the decay rate of Hi, th^ is the lifetime, and Ff^s = T^I^at/ + T^'^tv + 
F^i^TV + F^',^^, is the total absorption rate of Hi. 

Because the emitted high energy X° is not stopped in the thermal bath, its lifetime 
becomes longer by a factor of Ej^o /m^o due to the relativistic effect. Then, the decay rate 

is estimated by F^^^ = r~o m^^o /i?^o . Because the emitted -ft'^'s are distributed in the wide 
energy range, for convenience we computes the mean kinetic energy Ej^o which is obtained 



by weighting the kinetic energies for their distribution (see Fig. In Fig. |T0|, i^^^o is plotted 



as a function of 2£'jet. 



IV. BBN COMPUTATION IN HADRON-INJECTION SCENARIO AND 
COMPARISON WITH OBSERVATIONS 

In this section we perform the BBN computations in the hadron-injection scenario. Then 
we compare the theoretical prediction of the light element abundances with the observational 
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light element abundances. In the computations we assume that the massive particle X 
decays into three bodies (i^jet = tux/^) and two jets are produced at the parton level (i.e. 
the number of jets iVjet = 2). The above choice of a set of model parameters E-^ex, and iVjet is 
not unique in general and is obviously model dependent. For E-^ex, however, since we study 
the wide range of the mass, we can read off the results by rescaling the mass parameter. In 
addition, for the modification of A'jet since the second term in the right hand side in Eq. (|33D 
scales as oc iVjet, we only translate the obtained results according to the above scaling rule 
and push the responsibility off onto the number density nx- 

As we noted in the previous sections, it is a remarkable feature that the predicted "^He 
mass fraction Y tends to increase in the hadron injection scenario because ^He is the most 
sensitive to the freeze-out value of the neutron to proton ratio in the beginning of BBN. 
Since protons which are more abundant than neutrons are changed into neutrons through 
the strong interactions rapidly, the freeze out value of n/p increase extremely if once the net 
hadrons are emitted. In addition, D is also sensitive to the neutron number after T ^ O.lMeV 
because the free neutrons can not burn into ^He. 

To see the rough tendency, we plot the upper bounds for BhUx /s in Fig. ^ which come 
from each observational 2a upper bound for ^He, and D as a function of the lifetime tx 
at the baryon to photon ratio 77 = 5 x 10~^°. is the hadronic branching ratio of X, 
and nx/s denotes the number density of X per entropy density s. The mass is fixed to 
be typical value, mx = lOOGeV. From the figure, we find that for the shorter lifetime 
Tx ~ 10^^ sec, the hadron injections do not affect the freeze-out value of n/p and do not 
change any predictions of SBBN. However, if the lifetime becomes longer tx ^ 10~^ sec, the 
freeze-out value of n/p ratio is increased by the hadron- induced inter-converting interactions 
and the produced neutron increases the ^He abundance because most of the free neutrons 
burn into ^He through D. Then, nx/s is strongly constrained by the upper bound of the 
observational ^He abundance. For tx ^ 10^ sec, since the produced free D can no longer 
burn into ^He, the extra free neutrons still remain in D. Then nx/s is severely constrained 
by the upper bound of the observational D/H. For the constraint from high D, i.e. D/H 
< 3.0 X 10""^, we obtain the milder upper bound than low D because more productions are 
allowed from the observation. 

However, you can easily find that these constraints are obtained only when rj is fixed. If 
we chose the other r], e.g. which predicts more D/H than the upper bound of the observation 
in SBBN, then, the almost all parameter regions would have been excluded because both D 
and ^He tend to increase in the hadron-injection scenario. Namely, any constraints, which are 
obtained when we fix 77 a priori, have little meaning. To correctly compare each model in the 
various parameters {rj, tx, and nx/s), we should perform the Maximum Likelihood analysis 
and the fitting in wide parameter region including both the observational and theoretical 
errors. To estimate the theoretical errors, we perform the Monte Carlo simulation including 
the theoretical uncertainties which come from experimental errors of nuclear reaction and 



The Li abundance is mildly constrained from the observation and is much weaker than the 
others. In addition since it has a complicated dependence for r], we do not plot it here. Of course, 
however, we include ''Li in Monte Carlo simulation and Maximum Likelihood analysis which will 
be discussed below. 
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hadron-nucleon reaction rates. 

Concerning the detail of the executions, we have aheady explained in IIB . For the 
hadron-nucleon interaction rate, we adopt 50% error for each cross section because there are 
not any adequate experimental data for the uncertainties of cross sections. Therefore, we 
take the larger errors to get a conservative bound here. 

In Fig. |l2]we plot the results of the fitting at 95% C.L. in (rx, BhUx / s) plane projected 
on 7] axis in the case of low D which is obtained by performing the Maximum Likelihood 
analysis. The region below the line is allowed by the observations ("^He, D, and ^Li) for the 
various mass of X. If rrtx becomes heavier, more hadrons are emitted in the decay, and the 
upper bound becomes more stringent. Comparing the case of mx = 100 GeV in Fig. ^ with 



that in Fig. |Il|, the upper bound obtained in the Monte Carlo simulation is milder. That 
is because we did not adopt the naive 2a observational upper bounds with fixed 77, but we 
searched the wide range of rj, not forgetting tx and nx/ and we performed the Maximum 
Likelihood analysis to include both all the observational and theoretical uncertainties. In 



Fig. |T3| we also plot the results of high D. Compared to the case of low D (Fig. |T^), the 
obtained upper bound becomes milder because more D are allowed by the observations in 
high D case. 

As we also discussed in the previous section, the above treatment might underestimate 
the deuterium abundance for tx ~ 150 sec because deuterium is produced by the destruction 
of ^He by the high-energy free neutrons in such a relatively late epoch. Therefore, that means 
we obtained the conservative limits only for longer lifetime than 150 sec in this paper. 

Here we consider one of the concrete models of the hadronic decay. If we assume that the 
parent massive particle is gravitino and that it mainly decays into a photon and a photino 
(■^3/2 7 + 7)) the lifetime r3/2 is related to the gravitino mass m^i2 as 

r3/2^4xl0^secx (^^)"'. (36) 

In addition, if we assume that the gravitino is produced through the thermal scattering 
in the reheating process after inflationQwe relate the abundance 11^12/ s of the gravitino 
with the reheating temperature Tr 0, 



In Fig. [T^ we plot the upper bound on the reheating temperature after inflation at 95% C.L. 
as a function of the gravitino mass ^3/2- The solid line (dashed line) denotes the case of 
low D (High D). The region below the line is allowed by the observations. As we discussed 



^^In the last two years, it has been claimed that gravitinos might be also produced in the preheating 
epoch non-thermally pC|-32|. On the other hand, these days it was pointed out that such a effect is 
not important if we realistically consider two chiral multiplets to distinguish between inflatino and 
goldstino 1 33 1, although it may depend on the mixing of the supersymmetry breaking sector and 
the inflaton sector |34]. If the non-thermal production is effective, however, the obtained constraint 
might be severer. 
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before, Bh can become ~ 0{a) at least even if the main decay mode is only photons, because 
photon has the electromagnetic coupling with qq, i.e. {Bh = 0.01 - 1). For 171^/2 ~ lOTeV, 
they mean the conservative upper bound. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have discussed the effects of the late-time hadron injection on the 
primordial nucleosynthesis which are caused by the decay of an unstable massive particle X 
when the lifetime is relatively short 10~^ sec ^ tx ^ 10*^ sec. If the massive particle decays 
into quarks or gluons, they quickly fragment into hadrons. Then the high energy hadrons 
would be emitted into the electromagnetic thermal bath near the BBN epoch. Because the 
background photons and electrons are sufficiently energetic in the epoch, such high energy 
hadrons lose their almost all kinetic energies through the electromagnetic interaction, and 
they are approximately stopped before they interact with the background nucleons (p and n) 
except for neutral kaon K^. Then, they scatter off the background nucleons by the threshold 
cross sections only for the exothermic reactions and can extraordinarily inter-convert p and 
n each other strongly through the hadron-nucleon interaction even after the freeze-out time 
of the neutron to proton ratio n/p. At that time it is expected that the background proton 
tends to be changed into neutron through the strong interaction since protons are more 
abundant than neutrons, and n/p tends to increase. As a result, in particular, the abundance 
of ^He extraordinarily increases because it is the most sensitive to the freeze out value of 
n/p. Then, we can constrain the abundance of X and obtain the informations of tx from 
the observational light element abundances. 

Here we have studied the hadron injections by using the JETSET 7.4 Monte Carlo 
event generator |H| to quantitatively understand the hadron jets to agree with the collider 



experiments [^. Thanks to the treatments, we can estimate the number of the emitted 
hadrons as a function of the energy of jets, i.e. as a function of the mass of X, even 
in the regions where there is no experimental data. In addition we can also obtain the 
energy spectrum of the emitted for various masses of X. This is very important in the 
computations because is never stopped in the electromagnetic plasma, and we should 
know the energy distributions of -ft'^'s. On the other hand, we also have estimated the 
energy-dependent cross sections for iT^-nucleon scattering using the existing experimental 



data [27,29|. With these data, we could properly include the hadron-injection effects in BBN 
computations. 

To estimate the theoretical errors, we performed Monte Carlo simulation including the 
theoretical uncertainties which come from those of the elementary nuclear reaction rates 
and hadron-nucleon interaction rates. To obtain the degree of agreements between theory 
and observation, we performed the Maximum Likelihood method and the fitting including 
both the observational and theoretical errors. To correctly compare each model in the various 
parameters (rj, tx, and nx/s), the above procedure is quite crucial because a constraint 
which is obtained when we intentionally fix the parameters has little meaning. 

As a result, we have obtained the upper bound on the abundance nx/s as a function 
of the lifetime tx to agree with the observations for the wide range of the mass mx = 10 
GeV - 100 TeV which are relevant for various models of supergravity or superstring theory. 
However, we might have underestimated the deuterium abundances where the lifetime is 
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longer than C(IO^) sec because deuterium can be produced by the destruction of ^He by the 
high-energy free neutrons, i.e. "hadro-dissociation" effects which we ignored in this work. 
Therefore, if the process is effective, that would mean we obtained the conservative upper 
bounds only for tx ^ (9(10^) sec. In the separate paper, we will comprehensively study 
the subject [0. We have also applied the results obtained by a generic hadronic decaying 
particle to gravitino ip3/2- Then we have got the upper bound on the reheating temperature 
after primordial inflation as a function of the mass, Tr ^ 10'' — 10^ GeV (Tr < 10® — 10^ GeV) 
for ms/2 = 10 - 100 TeV at 95 % C.L. in the case of low D (high D). 
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FIGURES 
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1 10 100 
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FIG. 2. Plot of the averaged charged-particle multiplicity {Nch}- This represents the total 
number of the charged hadrons emitted per e~^e~ annihilation and per two hadron jets as a function 
of ^/s (= 2 -Ejet); where ^/s denotes the center of mass energy, and Ejet is the energy per one 
hadron jet. The solid line denotes the value obtained by using the JETSET 7.4 Monte Carlo 
event generator. The filled circle denotes the data points of e~^e~ collider experiments. Error is 
quadratically added for the statistical and systematic one. Here (A'ch) is defined as the value after 
both Ks and A*^ had completely finished to decay. 
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FIG. 3. Plot of the spectrum of the produced mesons (vr"'" + vr", + K^, and Kf^) as a 
function of the kinetic energy -Bkm- This is the case that the center of mass energy is y/s = 91.2 
GeV which corresponds to the resonance. They are computed by using the JETSET 7.4 Monte 
Carlo event generator. 
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FIG. 4. Plot of the spectrum of the produced baryons ((a) n + n and (b) p + p) as a function 
of the kinetic energy -Ekm- This is the case that the center of mass energy is ^/s = 91.2 GeV which 
corresponds to the resonance. They are computed by using the JETSET 7.4 Monte Carlo event 
generator. 
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FIG. 5. Plot of the averaged number of the produced hadrons as a function of 2£'jet(= \^), 
where -Ejet denotes the energy of one hadron jet. The number is defined by the value per two hadron 
jets. (A'ch) denotes the averaged charged-particle multiplicity (thick solid line). The number is 
obtained by summing up the energy distribution. The dotted line is 7r+ + vr^, the short dashed 
line is + K~, the thin solid line is K^, the dot-dashed line is p + p, and the long dashed line 
is n + n. They are computed by using the JETSET 7.4 Monte Carlo event generator. 
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FIG. 6. Plot of the distribution of K^^ produced in the e+e~ annihilation as a function of the 
kinetic energy. It is the case that the center of mass energy is ^/s = 91.2 GeV which corresponds 
to the resonance. It is computed by the JETSET 7.4 Monte Carlo event generator. 
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FIG. 7. Plot of the data of the cross sections as a function of the kinetic energy of the 
beam. 
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FIG. 8. Plot of the spectrum of the produced through the hadron fragmentation of qq pair 
emitted from e^e~ annihilation, x (= £'kin/\/s) denotes the normalized kinetic energy -Ekim and 
^/s denotes the center of mass energy of e^e~ collision. They are computed by using the JETSET 
7.4 Monte Carlo event generator. 



27 




FIG. 9. Plot of the averaged cross sections for p + — > p + • • • and p + -fC£ — > n + • • • as a 
function of the energy of two jets (= 2£'jet)- 
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FIG. 10. Plot of the mean kinetic energy of -fC£ which is obtained by weighting the kinetic 
energies for their distribution as a function of 2£^jet, where 2Ejet is the energy of two hadron jets. 
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FIG. 11. Plot of the rough upper bound of B^nx/s from the observational 2a upper bounds 
of *^He (solid line), and D (dashed line) for high D or low D as a function of the lifetime of the 
massive particle X. is the hadronic branching ratio of X, and nx/s denotes the number density 
of X per entropy density s. Here the baryon to photon ratio is = 5 x 10""^° and the mass of X 
(mx) is fixed to be 100 GeV. The observational upper bounds are obtained by adding the errors 
in quadrature. 
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FIG. 12. Plot of the contour of the confidenee level (C.L.) in [Tx^BhUx/ s) plane for low D. 
The region below the line is allowed by the observations at 95% C.L. tx is the lifetime of X, is 
the branching ratio into hadrons, and nx/s denotes the number density of X per entropy density. 
It is the case that the mass of X is, mx = 100 TeV (solid line), 10 TeV (dotted line), 1 TeV 
(dashed line), 100 GeV (long dashed line), or 10 GeV (dot-dashed line) respectively. 
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FIG. 13. Plot of the contour of the confidenee level (C.L.) in (tx , B^nx /s) plane for high D. 
The region below the line is allowed by the observations at 95% C.L. tx is the lifetime of X, is 
the branching ratio into hadrons, and nx/s denotes the number density of X per entropy density. 
It is the case of the mass mx = 100 TeV (solid line), 10 TeV (dotted line), 1 TeV (dashed line), 
100 GeV (long dashed line), or 10 GeV (dot-dashed line) respectively. 
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FIG. 14. Plot of the upper bound on the reheating temperature after inflation at 95% C.L. 
as a function of the gravitino mass 777.3/2- Here is the branching ratio into hadrons (= 0.01 - 
1). The solid line (dashed line) denotes the case of low D (high D). The region below the line is 
allowed by the observations. 
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